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Abstract. 

We estimate the relativistic binaries merger rate and redshift dis- 
tribution population synthesis of an ensemble of evolving close binaries. 
l/^ ' Results of such simulations definitely depend on the cosmic star forma- 

CN i tion rate history, which is different in galaxies of different types. This 

leads to a difference in merger rate in spiral and elliptical galaxies, which 
in principle can be an observational test for GRB models. Also a fit of 
C^ ' BATSE long bursts log A^-log P is performed, showing that only a wide 

lO . (2 — 3 orders of magnitude) luminosity function provides a good fit. 
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ON 

Q> ■ 1. Introduction 

^ ■ 
Q.. For many years the binary relativistic star merger remains one of the most 

valuable models of gamma ray bursts (Blinnikov et al, 1984, Paczyhski , 1991). 
Some recent alternative models, involving a collapse of massive rotating star 
^ I (Paczyhski , 1998), possibly easier fit the energy requirements for GRB, and are 

C^ . in nice agreement with the discoveries of GRB optical counterparts in actively 

star forming galaxies (Bloom et al, 1998). 

Nevertheless, the old merger model is not excluded, as the observed binary 
. , radiopulsars display the examples of definite merger precursors. 

C^ ■ This paper is devoted to the studies of close binary evolution which leads to 

the mergers of binary neutron star and black hole systems (NS+NS and NS+BH) 
which are suspected to be able to produce a GRB. 

The collapsar GRB and the merger GRB can have similar physical mech- 
anisms for energy extraction (Lee, 1999) and production of radiation (Piran, 
1999, Usov, 1999). Nevertheless, they can be discriminated by evolutionary 
considerations. 

2. The model of evolution 

We use the "Scenario Machine", initiated by Kornilov and Lipunov (1995), and 
later developed by Lipunov et al (1996), as a stellar evolution engine. The basic 
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evolution model is similar to one of Vanbeveren et al (1998), which are the best 
to reproduce the galactic population of massive binaries. 

During the evolution of a binary its orbital separation changes in various 
modes of mass and angular momentum transfer. The most dramatic change of 
the orbit takes place during a supernova explosion, or a collapse into a black 
hole. It is usually supposed that a newly formed compact object obtains a 
kick velocity due to some asymmetry in the collapse. To explain the observed 
velocities of radio pulsars, one should assume the kick velocity to be of the order 
of 200 - 600 km/s (Lyne, Lorimer, 1994). 

As the observations suffer from many selection effects, the exact shape of the 
distribution is not precisely known. So we perform the calculations for several 
values of the kick velocity, in the range from to 600 km/s, to show the influence 
of this parameter on the results. 

3. The event age distribution 

The life of the merging system consists of two important phases: first, the nuclear 
powered evolution of the normal stars, and second, the gravity wave powered 
orbit shrinking phase. The characteristic time scale of the first phase occupies 
the range from 3 ■ 10^ years for the most massive stars to ~ lO^years for the least 
massive stars being able to produce a neutron star. The gravitational inspiral 
time is determined by the parameters of the orbit after the formation of second 
compact object, and can be much greater than nuclear lifetime, being on average 
of the order of a billion years. 

The age of a merging binary is a sum of the nuclear lifetime and the grav- 
itational inspiral duration. The distribution of the ages of merging binaries is 
a direct output of the population synthesis procedure. Another sense of this 
distribution is the merger rate time history after a simultaneous ((5-like) star 
formation. It is shown in fig.|| for both NS+NS and NS+BH mergers. 

The merger age distribution shows a power-law behavior with a slope ~ 
—1. For NS+NS mergers, there is a strong dependency of the sharpness of the 
age distribution on the kick velocity. High kicks make the mergers on average 
younger. Without a kick, there are no mergers with ages less than 10* years. 

For NS+BH mergers, the presence of a kick velocity is very important: 
without a kick, all the binaries are very wide and have very long inspiral times, 
usually greater than the Hubble time. 

The power law distribution has a "heavy tail" : though a half of the mergers 
take place before several hundred million years after the star formation, there is 
also a significant part of mergers that take part after billions of years. 

4. The distribution of merger redshifts 

The distribution of merger redshifts can be obtained by a convolution of the star 
formation rate history and the merger age distribution. The present knowledge 
of the star formation rate (SFR) history is based on the works of Madau (1996). 
Earlier (J0rgensen et al, 1995) we assumed a simple two-parametric model of 
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Figure 1. The binary merger age distribution. Shown as the event 
rate after a simultaneous formation of a typical galaxy (IO^^Mq). Left: 
NS+NS, Right: NS+BH. 
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Figure 2. The merger redshift distributions (rate per unit z). 
NS+NS; Right: NS+BH. 
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SFR, containing an initial star formation burst (where e of all the stars were 
formed) at 2; ~ 5, and a subsequent constant SFR. 

In this work we base on the Madau observational SFR, but still introduce 
an initial star formation burst, which should be responsible for the formation of 
the elliptical galaxies and spheroidal components of the spiral galaxies. Then, 
according to Fukugita et al (1998), the value of e should be ~ 2/3. The adopted 
cosmological model was (B.q = 75 km/s/Mpc, 0^ = 0.3, Q.\ = 0.7). 

The obtained NS+NS and NS+BH merger redshift distributions are dis- 
played in fig. ^ for spiral and elliptical galaxies (Madau and burst-like SFR, 
respectively) and different values of mean kick velocity. 

Obviously, these event redshift distributions are different for spirals and 
ellipticals. This allows to propose a new observational test for the GRB pro- 
genitor. If we detect a GRB in an elliptical galaxy, it should be most likely a 
NS+NS and NS+BH merger and not a collapsar, because the population of the 
ellipticals is very old and does not contain massive normal stars. 
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Figure 3. The evolution of the ratio of merger rates in spiral and 
elUptical galaxies. Left: NS+NS; Right: NS+BH. 

Thus, future identifications of GRB host galaxies can solve the collap- 
sar/merger dilemma. 



5. The GRB brightness distribution — logA^-logP 

Another implication of the merger age distributions is an attempt to construct 
a GRB brightness distribution ( logA^-logP ). This is done by a classic pro- 
cedure (Weinberg, 1972) on the basis of the redshift distribution (fig. ^) and a 
luminosity function. Observations do not evidence directly for the spread of the 
GRB intrinsic luminosities, though some techniques allow to estimate it roughly 
(Petrosian, 1999). We assumed a power-law luminosity function characterized 
by the slope and range. 

For comparison with the observations, we have chosen the long (Tqq > 1.5s) 
and relatively bright (> 1 photon/cm^) bursts from 4th BATSE GRB catalog 
(Meegan et al, 1998). Realistic GRB spectra were taken to estimate the spectral 
K-correction (Oke, Sandage, 1968). 

The best fits are displayed in fig. I4 The log A^-log P distribution is shown 
in a differential form multiplied by P^ in order to outline the difference from 
the euclidean case, in which the distribution should look as a horizontal line. 

No acceptable fit can be obtained if luminosity spread is less than 2 orders 
of magnitude. In the best fit, the luminosity function spread is ~ 2.5 orders of 
magnitude, and the most distant observed GRB are at redshift ~ 4.5. 



6. Conclusions 

The computer simulations of the evolution of the possible GRB progenitors 
provide useful information which facilitates interpretation of the GRB statistics. 
In the present paper, we would like to outline two main conclusions: 

A new interesting method which can help to determine the nature of GRB 
progenitors is the determination of the ratio of the rates of bursts in spiral and 
elliptical hosts. For GRB produced by mergers, the rate in ellipticals should be 
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Figure 4. The best fits of the BATSE long bursts log iV-log P . Left: 
NS-FNS; Right: NS-^BH. Triangles: The best fit. Circles: Without a 
luminosity spread. Boxes: Without the spectral K-correction. Thick 
line represents the BATSE observations. 

only several times less than one in spirals, and increase with the redshift. For 
the massive collapsar GRB it is practically impossible to occur in an elliptical 
galaxy. 

The analysis of the observed BATSE log A^-log P distribution with taking 
into account the evolutionary effects shows that the the spread of GRB intrinsic 
luminosity function can not be less than two orders of magnitude. 
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